arXiv: 1503.01102v4 [cs.IT] 9 Mar 2016 


1 


Cooperative Base Station Coloring for Pair-wise 

Multi-Cell Coordination 

Jeonghun Park, Namyoon Lee, and Robert W. Heath Jr. 


Abstract —This paper proposes a method for designing base 
station (BS) clusters and cluster patterns for pair-wise BS 
coordination. The key idea is that each BS cluster is formed 
by using the 2nd-order Voronoi region, and the BS clusters are 
assigned to a specific cluster pattern by using edge-coloring for 
a graph drawn by Delaunay triangulation. The main advantage 
of the proposed method is that the BS selection conflict problem 
is prevented, while users are guaranteed to communicate with 
their two closest BSs in any irregular BS topology. With the 
proposed coordination method, analytical expressions for the rate 
distribution and the ergodic spectral efficiency are derived as 
a function of relevant system parameters in a fixed Irregular 
network model. In a random network model with a homogeneous 
Poisson point process, a lower bound on the ergodic spectral 
efficiency is characterized. Through system level simulations, the 
performance of the proposed method is compared with that 
of conventional coordination methods: dynamic clustering and 
static clustering. Our major finding is that, when users are dense 
enough in a network, the proposed method provides the same 
level of coordination benefit with dynamic clustering to edge 
users. 


I. Introduction 

A. Motivation 

Coordination among base stations (BSs) is a powerful ap¬ 
proach for mitigating inter-cell interference in cellular systems 
m . It has been shown that the sum spectral efficiency scales 
linearly with the signal-to-noise ratio (SNR) when all the BSs 
are coordinated m. In practice, however, coordination with 
a large number of BSs may not be feasible due to excessive 
overhead associated with the coordination, e.g., complexity, 
channel estimation and channel feedback E). One practical 
solution for implementing multicell coordination is to form a 
BS cluster so that a limited number of BSs are coordinated to 
control intra-cluster interference with a reasonable amount of 
overhead a, 0. Unfortunately, when clustering-based BSs 
coordination is applied in a static way 0, 0, or equivalently 
in a random way O independent to users’ conditions, the 
performance is mainly limited by unmanageable out-of-cluster 
interference a, a, 0. 

Dynamic clustering, where a BS cluster is dynamically 
made by users, was proposed in 0 to overcome the perfor¬ 
mance limitation of static clustering. While dynamic clustering 
improves the signal-to-interference and noise ratio (SINR) of 
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a user significantly, there are several challenges that make it 
difficult to realize. In particular, different users who share the 
same spectrum can select the same BS at the same time. We 
refer to this as “the BS selection conflict problem.” To resolve 
this, users in a network need to send feedback regarding 
desirable BS clusters and a global scheduler is required to 
schedule BSs and users based on this feedback information. 
For this reason, performance chiefly depends on the choice of 
this scheduler. The scheduling itself, however, requires a lot of 
feedback overhead from the users and a substantial amount of 
coordination of all the BSs and the users, which causes huge 
system overhead. 

To solve the BS selection conflict problem without requiring 
global coordination, semi-static clustering was proposed in 0, 
Eol. In this method, multiple predefined BS clusters are used 
with different time-frequency resources. A set of BS clusters 
that use the same time-frequency resource is referred as “a BS 
cluster pattern.” For instance, in a square grid BS topology, 
four BS cluster patterns are required to cover the whole region 
without the BS selection conflict problem. Using the semi¬ 
static clustering, users can communicate with their favorable 
BS set. The main limitation of the existing work 0, cni is 
that the BS clustering methods are applicable for a regular 
BS topology where each BS is located regularly on a grid. 
In practice, however, the BS topology of cellular networks is 
irregular mi. In such case, BS clustering pattern design is 
challenging. In this work, by leveraging geometry properties, 
we make progress in designing BS clusters and BS cluster 
patterns in a semi-static way, which is applicable for any 
wireless networks whose topology is irregular. 

B. Related Work 

Cooperative strategies for dynamic clustering have been 
extensively studied in the literature 0, lfT^ - ll2^ . In 0, ifT^ - 
M, an algorithm for making a BS cluster in dynamic clus¬ 
tering was proposed. For instance, given the users’ conditions 
such as distances or channel conditions, an algorithm selects 
BSs to include into a BS cluster so that the joint achievable 
throughput of the formed BS cluster is maximized. In ifTSl - 
El, a resource allocation method to enhance sum spectral 
efficiency was proposed for dynamic clustering. By jointly 
allocating resources to each BS cluster depending on their 
channel qualities, the sum spectral efficiency in a BS cluster 
can be improved. In lfT8l - ll20l . a beamforming algorithm was 
proposed in dynamic clustering aiming complexity reduction 
or sum spectral efficiency improvement. In this work 0, 
mi-iiioi, a toy network, where only limited number of BSs 
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and users are assumed, was considered. In ||2TI| - Il24l . the 
performance of dynamic clustering was characterized in a 
random network where BSs are distributed according to a 
homogeneous Poisson point process (PPP). In 123, pair-wise 
dynamic clustering with rate splitting method was proposed 
assuming a Poisson network. In 123, dynamic clustering in 
two-tier network was proposed and analyzed also in a Poisson 
network. The main limitation of the prior work considering 
BS clustering in a random network m-M is that they did 
not provide a solution for resolving the BS selection conflict 
problem, rather they focused on the performance characteriza¬ 
tion (e.g., the coverage probability) under the assumption that 
a global scheduler allows such BS clustering to a typical user 
without the BS selection conflict problem. Global scheduling 
has higher complexity as overhead and scheduling are required 
between all possible BS clusters. The main difference between 
this paper and the prior work in lISTl - ll^ is that we focus on 
providing a systematic solution for the BS selection conflict 
problem without using a global scheduler, so that the proposed 
method is able to be implemented in a network with reasonable 
overhead. 


C. Contributions 

In this paper, we propose a method for designing BS clusters 
and BS cluster patterns for pair-wise multicell coordination 
in a cellular downlink network whose BSs’ locations are 
irregular. The goal of the proposed method is to ensure 
that selected users communicate with their two closest BSs 
avoiding the BS selection conflict problem. To do this, we 
leverage the concept of 2nd-order Voronoi region. Suppose two 
BSs are located at and dj € in a two dimensional 

space. Then, the set of closer points to d^ and dj than any 
other points is defined by the 2nd-order Voronoi region as 

V2 (di,dj) 

d e 

The core feature of the 2nd-order Voronoi region is that 
it characterizes an area where any point in that area has 
the minimum distances to {d,, dj}. If a BS pair located at 
{dijdj} forms a BS cluster and serves users in V 2 (di,dj), 
those users are guaranteed to communicate with their two 
closest BSs. Applying this to the whole network, the first step 
of the method is tessellating a network plane into the 2nd-order 
Voronoi regions of every feasible pair of BS locations and 
making a BS cluster according to 2nd-order Voronoi regions. 
Next step is designing BS cluster patterns by allocating time- 
frequency resources to each BS cluster for preventing the BS 
selection conflict problem. To do this, we adopt the concept of 
Delaunay triangulation and edge-coloring in graph theory. At 
the best of authors’ knowledge, it is the first work to solve the 
BS selection conflict problem in irregular BS topologies. By 
solving this, users are able to obtain BS coordination benefit 
by communicating with their two closest BSs. We explain the 
proposed idea more rigorously in Section II. 
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With the proposed method, we characterize the performance 
of the proposed BS coordination strategy in a fixed irregular 
network model. Under the premise that multi-user coordinated 
beamforming (CBF) EU is applied to mitigate intra-cluster 
interference, we derive analytical expressions for the rate 
coverage probability and the ergodic spectral efficiency as a 
function of relevant system parameters: 1) distances between 
BSs and users, 2) the number of antennas per BS, 3) the 
number of selected users, 4) signal-to-noise ratio (SNR), 
and 5) the pathloss exponent. Next, considering a random 
network model where a location of each BS is modeled by 
using a homogeneous PPP, we derive a lower bound on the 
ergodic spectral efficiency as a function of the relevant system 
parameters. 

To demonstrate efficiency of the proposed method, we 
use a system level simulation to compare the edge users 
sum throughput of the proposed method and that of other 
conventional coordination methods: dynamic clustering and 
static clustering. One main finding is that when users are dense 
enough, the proposed method achieves the same edge user 
performance as dynamic clustering in which a global scheduler 
is assumed. This gain comes from the fact that the proposed 
method ensures that selected users communicate with their two 
closest BSs with high probability as the user density increases. 

The remainder of the paper is organized as follows. In Sec¬ 
tion II, the system model including the proposed cluster pattern 
is explained. In Section III, the rate coverage probability is 
characterized and the ergodic spectral efficiency is obtained in 
Section IV. In Section V, the ergodic spectral efficiency in a 
random network is analyzed, followed by Section VI which 
provides comparisons with simulations. Section VII concludes 
the paper. 


II. System Model 

In this section, we introduce the system model used in 
this paper. We first describe the network model, the pair¬ 
wise coordination model, and user association assumptions. 
Then the proposed BS cluster and BS cluster pattern are 
explained by using a toy example, and extend them to a general 
network model in the next subsection. Further, we introduce an 
algorithm for enhancing the proposed method, called as “edge 
cutting algorithm.” Signal model and performance metrics are 
specified in the following subsection. 

A. Network Model 

We consider a fixed downlink cellular network where each 
BS has N antennas. As illustrated in Fig.[T](a), BSs are located 
at irregular points on a two dimensional plane. We denote 
each location of a BS as di G for i € so that M = 
{di| i G Z+} includes all the locations of the BS on the plane. 
A BS located at d^ for i G has its own Ist-order Voronoi 
region Vi (d^), and each Ist-order Voronoi region tessellates 
the network plane into |A/"| regions. 

B. Pair-wise Multicell Coordination 

In this paper, we focus on the case where only two adjacent 
BSs form a cooperative cluster. The rationale behind focusing 
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(a) (b) 

Fig. 1. (a) An example of the considered network model. This example 

is generated by a repulsive point process where a distance between any two 
BSs is larger than 1.4km. (b) A graph induced by the Delaunay triangulation 
corresponding to (a). 

on pair-wise BS coordination is that it is highly representative 
for characterizing coordination gain without causing too much 
overhead. For example, in ETll . it was shown that a BS 
cluster including more than two BSs decreased the ergodic 
spectral efficiency when considering signaling overhead. This 
is because including more BSs into a BS cluster requires more 
overhead for estimating channel coefficients of the BSs in the 
cluster, which decrease the ratio of the transmitted data in 
a packet. This degrades the spectral efficiency. For pair-wise 
multicell coordination, a BS located at dj is able to coordinate 
with a BS located at dj if A (V 2 (d^, dj)) ^ 0, where A (■) is 
the Lebesgue measure in two dimensional space. 

C. User Association 

We assume that 2K single antenna users are selected in each 
2nd-order Voronoi region irrespective of the size of the region. 
Among them, K users are associated with the BS composing 
the corresponding 2nd-order Voronoi region while the other 
K users are associated with the other BS. For instance, if 2K 
users are in V 2 (d^, dj), K users in Vi (di) n V 2 (di, dj) are 
associated with the BS located at di and the other K users in 
Vi (dj)nV 2 (di, dj) are associated with the BS located at dj. 
This underlying assumption is justified by assuming that the 
number of users in a network is far larger than the number of 
deployed BSs, so that there are at least 2K users in each 2nd- 
order Voronoi region satisfying our assumption. More specific 
user scheduling algorithms are beyond the scope of this paper. 

D. The Proposed Clustering Model - Toy Example 

In this subsection, we explain the proposed clustering model 
with a motivational example. To this end, suppose that BS i 
is located at d^ in a finite area .4 C for i G {0,1, 2, 3}, 
as illustrated in Fig. |2] By exploiting the notion of the 2nd- 
order Voronoi region O, it is possible to create ( 2 ) = 6 
different pair-wise cooperative areas V 2 (di,dj) for i ^ j and 

i, j G {0,1,2,3} such that all points within the given region 
V 2 (di,dj) have the same two nearest BSs, i.e., BS i and BS 

j. In the proposed algorithm, users located in V 2 (di,dj) are 
served by a BS pair located at {d^, dj}, so that they are guar¬ 
anteed to communicate with the two closest BSs. Next, once 
each of 2nd-order Voronoi region is defined, multiple cluster 
patterns Vi for i = {1,..., Lj are created, each of which con¬ 
tains different 2nd-order Voronoi regions and 
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Fig. 2. A motivational example where all BSs are assumed to be on a finite 
triangular plane A, each of which is located at for i G {0,1, 2, 3}. The 
whole plane can be divided by the 2nd-order Voronoi region of {d^, d^} for 
i,j G {0,1, 2, 3} and i 7 ^ j. The areas which has the same shade pattern are 
assigned to the same cluster pattern, resulting in that three cluster patterns, 
i.e., 7 ^ 1 , 7^2 5 and Vs, ai'e defined. 

Each cluster pattern uses different time-frequency resources, 
so that there is no BS conflict between two different cluster 
patterns. An example of the BS selection conflict problem 
occurs when V2(do,di) G Vi and also V2(do,d2) G Vi. 
In this case, BS 0 is conflicted, and it cannot serve all the 
associated users. To avoid this, one possible example is to 
create six cluster patterns such that Vi = {V 2 (do, di)}, 
V 2 = {V2(do,d2)}, V 3 = {V2(do,d3)}, Vi = {V2(di,d2)}, 
V 5 = {V 2 (di,d 3 )}, and Ve = {V 2 (d 2 ,d 3 )}. This example, 
however, requires too many time-frequency resources, i.e., 6. It 
is even larger than the case where each BS uses different time- 
frequency resources without BS coordination, i.e., 4. A better 
example is to create three BS cluster patterns such that Vi = 
{V2(do, di), V2(d2, ds)}, V 2 = {V2(do, d2), V2(di, d3)}, 
and V 3 = {V 2 (do, d 3 ), V 2 (di, d 2 )}. With this example, BS 
conflicts are prevented and the required number of time- 
frequency resources is reduced to 3. To figure out how much 
performance gain is achieved by using this cluster pattern, 
assume that a user is in cell-edge region where the interference 
chiefly comes from the closest interfering BS. When no BS 
coordination is applied, the performance of the user might 
be severely degraded by the dominant interference which has 
almost same power with the desired signal. When each BS 
has its own separate resource, the dominant interference can 
be mitigated but it requires 4 time-frequency resources. By 
applying the proposed pattern, the dominant interference is 
mitigated with 3 time-frequency resources, which can lead to 
the performance improvement. 

E. The Proposed Clustering Model - General Network 

We extend the proposed clustering concept to a general 
network. At first, we illustrate in Fig. [3] how the proposed BS 
clusters are formed in a general network. Focusing on a BS 
located at do and its adjacent BSs at dj for j G {1, 2,..., 5}, 
users in V2(do,dj) for j G {1,2,..., 5} are served through 
a BS pair (dp, dj} by using the proposed clustering method. 
After forming BS clusters, we design cluster patterns V^ for 
£ G {1,2,...,L} where V^ = { V2 (di,dj)| i, j G Z+j. From 
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Fig. 3. An illustration of cooperative areas characterized by the 2nd-order 
Voronoi region (denoted as dotted red line) focusing on the BS located at 
do. Each cooperative area marked by different shade patterns is assigned to 
a different cluster pattern to avoid BS conflicts. 

the intuition of the toy example, we summarize the conditions 
for a desirable cluster pattern. 

1. Vi = (all users in the network are covered). 

2. Any two V 2 (dj,dj) e 'P^, 122 (d^, d^,) S Vt for ^ e 

L}, i ^ j ^ V ^ w. (to avoid BS conflicts). 

3. Small L satisfying 1) and 2): (to minimize the use of 
unnecessary resources). 

In a network whose BSs’ locations are irregular, it is chal¬ 
lenging to figure out how to map each of the BS clusters into 
appropriate BS cluster patterns satisfying the above conditions. 
To solve this, we draw a graph G (Af) whose a vertex is 
di G Af and an edge is made by Delaunay triangulation 
defined in JV. Delaunay triangulation for Af is a triangulation 
of a plane, where each vertex of a triangle is d^ G M. 
One important condition of the Delaunay triangulation is that 
when drawing the circumcircle of a triangle, there should be 
no point of M inside that circumcircle. Note that Delaunay 
triangulation always can be defined if there are more than three 
points, i.e., |Af| > 3. The example of G (A/") is illustrated in 
Fig. [U (b), when M is corresponding to Fig. [T] (a). A key 
relationship between 2nd-order Voronoi region and Delaunay 
triangulation is that if there is an edge between d^ and dj, 
i.e., E{di,dj) = 1 in G{M), the corresponding 2nd-order 
Voronoi region is non-empty, i.e., A(V 2 (di,dj)) ^ 0 ll28ll . 
For this reason, designing BS cluster patterns, i.e., mapping 
V 2 (di,dj) into a cluster pattern, is equivalent to mapping 
E (di, dj) into a cluster pattern. When considering the cluster 
patterns as “colors,” the cluster pattern design problem is 
equivalent to the edge-coloring in graph theory 1^ . The goal 
of edge-coloring is to assign colors to the edges of the graph 
so that any two edges sharing the same vertex have different 
colors, while minimizing the required number of colors to 
cover the whole graph. By using edge-coloring for G (Af), 
we are able to create cluster patterns Vi, £ G {1, 2,..., L} that 
satisfies the above conditions. 

Now we study the relationship between the cluster pattern 
design and the edge-coloring specifically. At first, we char¬ 
acterize L, which is the required number of different time- 
frequency resources to avoid the BS selection conflict problem. 
Thanks to the equivalence between designing cluster patterns 
and the edge-coloring, L is equal to the minimum required 



Fig. 4. Illustration of two different Delaunay triangulations. The membership 
of the cluster for its cluster pattern is indicated by the edge label, (a) An 
example of the Delaunay triangulation of a symmetric network, i.e., the 
maximum degree is equal to the minimum degree, (b) An example of the 
Delaunay triangulation of an asymmetric network, i.e., the maximum degree 
is strictly larger than the minimum degree. 

number of colors for the edges. The minimum required number 
of colors is characterized in the following Theorem. 

Theorem 1 (Vizing’s theorem, |I30|). A simple planar graph 
G of maximum degree A(G) has chromatic index A(G) or 
A(G) -f 1 in general. 

Proof: See the reference IMIl . ■ 

The chromatic index is known as the required minimum 
number of colors and the degree means the number of edges 
connected to the corresponding vertex. In general, it is NP- 
complete to find the chromatic index given an arbitrary graph, 
but we can observe that a graph G made by Delaunay 
triangulation has chromatic index A(G) in many cases. Now 
we provide two examples to show how to create the cluster 
patterns for a given graph. The first example, described in 
Fig. m (a), is a symmetric network where every vertex of 
G (Af) has the same degree, 6. This example is fitted to a 
scenario where each BS is deployed with a guard distance 
for preventing that two BSs are located very close together. 
Specifically, each BS is randomly located within a certain 
circle, and a center of each circle is located at a grid. For 
this example, L = 6. According to the description in Fig. 0] 
(a), the cluster pattern Vi for f G {1,2,..., 6} can be designed 
as 


Vi -- 

= {V 2 

(do, di), V2 

(d 2 , dio), V2 (do, d 4 ), • 


V 2 -- 

= {V 2 

(do, d2), V2 

(do, dio), V2 (d 4 , dll), 

...} 

Vz-- 

= {V 2 

(do, do), 122 

(di, d 2 ), V2 (d 4 , do), • • 


V 4 -- 

= {V 2 

(do, d 4 ), 122 

(d 2 , do), V2 (do, do), • • 

•f 

V 5 -- 

= {V 2 

(do, ds), V2 

(di, do), V2 (d 7 , dg), • • 


Ve-- 

= {V 2 

(do, do), V2 

(d 4 , dg), V2 (do, dll), • 



by using edge-coloring. By allocating different time-frequency 
resources to each cluster pattern Vi for £ = {1,2, ...,6}, no 
BS conflict occurs and users are guaranteed to communicate 
with their two closest BSs set. 

Now we look an example in Fig. 0] (b), which describes 
an asymmetric network, i.e., each vertex of G (Af) can have 
different numbers of degrees. This example is fitted to a case 
where each BS is deployed with a more random manner than 
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the symmetric network case. In this case, L = 7 is enough. 
Similar to the previous example, the cluster pattern Ve for 


{1, 2,..., 7} can be 



iPl ={V2(do,d2), 

V 2 (ds, dg), V 2 (do, dj), • 


V 2 = {V2(d2,d3), 

V 2 (do, d 4 ), V 2 (do, dg), • 


'Ps = {V 2 (do, ds), 



iP4 = {V2(d4,d5), 

V 2 (do, dg) ,•••}, 


7^5 = {V 2 (d7, dg), 



iP6 = {V2(d2,d4), 

•••}, 


P? = {V 2 (da, d4), 




also by using edge-coloring. No BS conflict occurs as long as 
a different time-frequency resource is assigned to a different 
cluster pattern. In this paper, we do not specify an algorithm 
for solving edge-coloring, rather we assume that edge-coloring 
is solved perfectly. Some useful algorithms for edge-coloring 
are found in OTll . Il32ll . In OTIl . an edge-coloring algorithm for 
simple graphs was proposed and it was proven that the algo¬ 
rithm complexity is 0{A{G) lA/"! \/\JV\ logdAfl)). In ll^ . an 
edge-coloring algorithm for planar graphs whose A(G') > 9 
was proposed and the algorithm complexity was proven to be 

om v/|Ar|iog(|Ar|)). 

Observe that although the general edge-coloring problem 
is NP-complete, complexity is not a bottleneck since the 
proposed method (making BS clusters and BS cluster patterns) 
does not have to be performed frequently. Once BS clusters 
and cluster patterns are made, they are preserved until the ge¬ 
ometry of the BSs changes. Conventionally, the BS geometry 
would not be changed on the order of months or years. 

F. Edge Cutting 

In some BS topologies, particularly in an asymmetric net¬ 
work case, it is possible that unnecessary time-frequency 
resources are wasted. For instance, in Fig. m (b), we found 
that L = 7 time-frequency resources are required for the 
corresponding edge-coloring. Focusing on BS located at do 
in Fig. m (b), however, it only uses 5 resources therefore 
2 resources are wasted. Due to these wasted resources, the 
network sum throughput even can be degraded. To relieve 
this, the edge cutting algorithm is proposed. The key idea of 
this algorithm is decreasing the maximum degrees A(G) to 
Aec(G) < A(G) by cutting some edges in the original graph 
G(A/’). Specifically, in G (Af), if a vertex G A/” has degree 
D (di) > Aec(G), D (di) — Aec(G’) edges are selected and 
cut. Repeating this to every vertex in a graph, we construct a 
cut-graph G'^“* (Af) whose the maximum degree is Aec(G')- 
Note that the constructed cut-graph G'^“* (Af) is a subgraph of 
G (A/"). Assuming that an edge E (d^, dj) is cut in G'^“* (A/"), 
users in V 2 (d^, dj) are not served. By doing this, the required 
time-frequency resources are Aec(G). By applying the edge¬ 
cutting algorithm, the overall network performance can be 
improved by saving time-frequency resources. Specifically, the 
required time-frequency resources L decreases to Aec(G) 
from A(G) (See Theorem [Tli, where Aec(G) can be con¬ 
trolled as an algorithm parameter. 


There are two noticeable points in our edge cutting algo¬ 
rithm. The first one is that selecting the edges to cut has 
a significant effect to the performance of the algorithm. A 
desirable way is selecting edges inversely proportional to 
the area of the corresponding 2nd-order Voronoi region. For 
instance, assume that we have two candidate edges for cutting, 
E{di,dj) and £'(di,dfc). If A(V 2 (di,dj)) > A(V 2 (di,dfc)), 
we cut i?(di,dfc) since its 2nd-order Voronoi area is smaller 
than the other’s area. By doing this, a large cooperative area 
remains while saving the resources. The second point is that 
we can restore edges after finishing the edge cutting. The 
intuition is that the ideal situation for the proposed method is 
that every BS has the same degree, as in a symmetric network 
case Fig. 4 (a). For this reason, we find a vertex whose degree 
is less than Aec(G) and re-connect edges while satisfying the 
maximum degree of the graph is Aec(G). The edge restore 
order is reverse of the edge cut order, so that E (di,dj) is 
restored if A(V 2 (d^, dj)) > A(V 2 (d^, d^)). The pseudo code 
description of the edge-cutting algorithm is described in the 
following using a MATLAB style notation. 


Algorithm 1 Edge-cutting algorithm 
Input: G(Af), Nbs = |Af|, Aec(G). 

Initialize: Let S'Neighbor(di), * G A/" be a set of vertices 
connected via the edge to d^ in G(A/") 

(Area-based edge cutting) 
for i = 1 : Nbs do 

if |*5’Neighbor(di) I > AEc(f^) thCU 
Find dj G 5’Neighbor(dj) 

Cut dj in order of (A(V 2 (d,, dj))) (ascending order) 
so that |5'Neighbor(di) I ^ Aec(G) 

end if 
end for 

(Area-based edge restore) 
for i = 1 : Nbs do 
if |»5’Neighbor(di) I ^ Aec(G) then 

Find dj cut in the area-based edge cutting 
Restore dj in order of (A(V 2 (d^,dj))) (descending 
order) satisfying max(5'Neighbor(dfc)) < AEc(G),fc G Af 
end if 
end for 

Output: G‘='^* (AC) 


Remark I. Since users in the cut region are not served in 
the proposed algorithm, there can be a fairness issue. This 
is able to be resolved by using the proposed method and 
conventional single cell operation method alternatively. For 
example, we can use the proposed method during time fraction 
a and conventional single cell operation during time fraction 
(1 — a). Every region is covered during time fraction (1 — a). 
In addition to that, more advanced scheduling algorithm such 
as proportional fairness can be also incorporated. 

Remark 2. In the proposed edge cutting algorithm, Aec(G) 
can be determined in various ways. One promising way is to 
make that all vertices in G{N') have the same degree, so that 
time-frequency resources are equally allocated to every BS. 
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To do this, Aec(G') should be equal to the minimum degree 
of G, i.e, minZ)(di), Vd^ € A/". Other ways to determine 
Aec(G') can be incorporated depending on BS geometrical 
features and user density. For instance, one can formulate an 
optimization problem as in /Ea?, where Aec(G) is maximized 
while guaranteeing the difference between maxD(di), Vd^ € 
jV and minZ?(di), Vd^ S M is less than a particular value. 

Finally, we summarize the procedure of the proposed 
method including the edge cutting algorithm as in the fol¬ 
lowing; 

1. Tessellate a network plane into 2nd-order Voronoi regions 
V 2 (di,dj) for dj,dj e J\f. 

2. Make a BS cluster {di, dj} according to the corresponding 
V 2 (di, dj) for dj, d^ G M. 

3. Draw a graph G(J\f) by Delaunay triangulation. 

4. Perform the edge cutting algorithm. 

5. By using edge-coloring, make BS cluster pattens Vi for 
i G {1,.., L} and assign different time-frequency resource 
to each BS cluster pattern. 

Now we investigate computational complexity of each step. 
Assuming we have lA/”! points, constructing \ff\ points 2nd- 
order Voronoi region requires 0(4 lA/”! log lA/”!) Il34l . drawing 
Delaunay triangulation demands OdA/”! log lA/”!) 135], and the 
edge cutting algorithm needs OdA/”! A(G)). In general, edge¬ 
coloring is NP-complete. The complexity is not important part 
in the proposed method since BS clusters and cluster patterns 
made by the proposed method are preserved until the geometry 
of the BSs changes, which occurs on the order of months or 
years conventionally. 


G. Signal Model 


Now we describe the signal model. We focus on a user 
associated to the BS located at do. We denote this user as the 
tagged user. We also assume that the tagged user is located on 
the origin. This assumption is easily generalized to arbitrary 
user location by shifting each BS’s location to — u for 
i G Z+, where u is a location of the tagged user. Under 
the assumption that the tagged user is in V 2 (do,di) and 
V 2 (do, di) is not cut in the edge cutting algorithm, a BS pair 
{do,di} forms a cluster according to the proposed cluster 
pattern. Let’s also assume V 2 (do,di) G Vi. In each BS 
located at d^, an information symbol vector sf G where 
sf = [sf sf jf] for i G is encoded to be transmitted 
to K users. The average power of this transmit symbol vector 


satisfies E 


< P. When transmitting the symbol vector. 


where h.[ G is the channel coefficient vector from the 
BS at di to the tagged user. The channel coefficients are 
assumed to have independent and identically distributed (IID) 
complex Gaussian entries with zero mean and unit variance, 
i.e., CA/^(0,1). The distance dependent pathloss is with a 
reference at Im and has an exponent of /3 and the additive 
Gaussian noise is rf that follows CM (O, cr^). A/f C A/" is a 
set of BSs included in the same cluster pattern with the BS at 
do, i.e., ffi = {dj;| V 2 (di, dj) G Vi,\/j}. The received signal 
can be separated into three parts: the desired signal and the 
inter user interference from the associated BS at dg, the intra¬ 
cluster interference from a BS set di, and the unmanageable 
out-of-cluster interference, respectively. 

Before transmission of an information symbol, pilot sym¬ 
bols are exploited to learn the downlink channel coefficient 
vector ho and h^. Then, the user sends the obtained CSI back 
to the corresponding BSs via an error-free feedback link. CSI 
at the transmitter (CSIT) within the same BS cluster is allowed. 
To mitigate intra-cluster interference, CBF is employed by 
using obtained CSIT. Assuming that our tagged user’s index is 
k, the beamforming matrix ’Vq and are designed to satisfy 
the following condition. 


maximize : 


subject to : 


(hO^vf 


(hO' 


0 ,k 


•' 0 ,k' 


'l.k'' 


= 0 for k' k 
= 0 for 1 < k" < K, 


(5) 

(6) 
(7) 


where the constraints of multi-user CBF are for nullifying 
the inter-user interference and for removing intra-cluster in¬ 
terference, respectively. Under the conventional zero forcing 
(ZF) conditions, it is always possible to find such Vq and 
Vj if 2K < N with high enough probability if each channel 
coefficient is IID. 


//. Performance Metrics 

If perfect CSIT within the same BS cluster is allowed, then 
the inter-user interference and intra-cluster interference terms 
in (IHi are zero. Assuming that indicates a modified channel 
coefficient after multiplying with a beamforming matrix Vq 
and Sq q is the information symbol for the tagged user, we 
have the following modified received signal 

dj€Afe\{dQ,di} 


a linear beamforming matrix Vf = [vf 


l! ■ 


U.k 


G and 
observation from 

l|do 


U,k 


G,K\ 


whose 


( 8 ) 


— 1 for k — 1,...,K is used. The The SINR of the information symbol s\ 


= 


the tagged user is given by 

(h^Q)^ V^s^o +\\di\\-f^/^hifvW 


1^1 SINK,. = 


0,0 


IS 


desired signal + inter user interference intra cluster interference 


Sd, 


eVA{do,di} 


VayPlK 


+ E (4) 

djeA/i\{do,di} 


out of cluster interference 


2 


riQ 

/^+ ||dof/SNR/iT’ 


(9) 
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(h')’'V; 


where It = Ed,€W,\{do,di>(il‘*)ll/ll‘*oll) 
and Pja^ = SNR. Given the system assumptions, the rate 
coverage probability is given as 

(SNR, lldoll, Vt, N, K, /3, 7 ) = P [log 2 (l + SINR,,) > 7 ] 

( 10 ) 

where Vg = {| G A/£\{do, di}} and 7 is the rate 
threshold. The ergodic spectral efficiency is then 

(SNR, lldoll, Vi, N, K, L, (3) = [log 2 (l + SINR,,)] . 

( 11 ) 

The pre-log term 1/L is used because L time-frequency 
resources are used in a network. L is determined as an 
algorithm parameter in the edge cutting algorithm. 


III. Rate Coverage Analysis 

In this section, we derive the rate coverage probability in 
a closed form under the interference limited assumption. To 
this end, we first introduce the following lemma. Lemma [T] 
provides distributions of the desired channel gain and the 
interference power. 

2 


the Laplace transform, i.e., E = (—1)™ ^ 

(O is represented as 


P/, 


N—2K 1 \m 


m—0 
N-2K 


m\ 


^ ^ (2-^ - 1)^ ^ {s) 


m—0 


m\ 


ds^‘ 


(14) 


S=2T-1 


Lemma 1. The desired channel gain 


where (s) is the Laplace transform of R, defined as 

^l(s)= TT I - - - 7 ) > ( 15 ) 

which completes the proof. ■ 

Computing the rate coverage probability ([I2]i is not 
easy since computing the derivative of the product of 

^1/1 -f s (||dj|| / ||do||)~^^ results in many terms. To pro¬ 
vide a simpler form of the rate coverage probability, the 
following Corollary is used to approximate the rate coverage 
probability. 

Corollary 1. The rate coverage probability (112b can be 
approximated as in the following form. 


follows Chi- 


squared distribution with 2 {N — 2K 1) degrees of freedom. 
The out-of-cluster interference R follows weighted sum of Chi- 
squared distribution with 2K degrees of freedom. 

Proof: See ET] and reference therein. ■ 

By using this Lemma, we derive the rate coverage. The 
following Theorem is the main result of this section. 

Theorem 2. In the interference limited regime, when the BSs 
form BS clusters by using the proposed method, the rate 
coverage probability of the tagged user is 

P‘ (SNR 00 , lldoll, Ve, N, K, /3, 7 ) = 

d' 


(SNR ^ 00, lldoll, Vi, N, K, P, 7) 

N—2K - -.m ^ ~ ^ 

i: siLiP (- 1 )" 


m—0 


ml 


ds^‘ 


djei>^\dmin I l|do|| 


-0 ^ 


K 


V 


/ lldmlnll V 

lldoll ) 


(16) 


where s = 2 '’' — 1 and dmin = argmin ||dj| 

dj 


N—2K - N m 


m—0 


ds^ 


n 

dj ^T)g, 


/ 


K 


( 1 \ 


1 


n 

I 

7 + s 

1 

0 

’ 

7 »(s)'v 


Proof: Using the Tailor expansion of the exponential 
don = l + ^ + §r + '’'^ when x is small enough, 
! 1 -I- X, therefore e~^ « 1/ (1 -|- x). Then we have 

K 


( 12 ) 


where s = 2 '*' — 1 and 7 is the rate threshold. 

Proof: As defined in (fTOb . the rate coverage probability 
is rewritten as 

P [log2 (1 + SINR|^) > 7] 

(“) , 


R 


> 12^ - 1 ) 



(b) 


( 2 ^- - 1 )" 

m—0 


ml 


{hr 




(13) 


where (a) comes from that ||d||^/SNR 0 when SNR 

00 , and (b) follows Lemma [T] and the complement cumulative 
distribution function of a Chi-squared random variable with 
degrees of freedom 2 {N — 2K -\- 1). Using the derivation of 


_i_ / ||dmin II \ 

- + lldoll ) 

which completes the proof. ■ 

Using the Taylor expansion, we approximate a product of 

( 1/1 -I- s (||dj|| / ||do||)~^j in an exponential form, so that 
we are able to calculate tne derivative of (fTbl) more simply 
than that of (fTTI) . 

Now we verify the approximated form of the rate coverage 
probability (fTbl) by comparing to the Monte-Carlo simulation 
results. For the verification, we use a network model described 









































- N=4,K=1 (simulation) 

□ N=4,K=1(analytical approximation) 

-N=3,K=1 (simulation) 

^ N=3,K=1 (analytical approximation) 

-N=4,K=2(simulation) 

^ N=4,K=2(analytical approximation) 

-N=10,K=5(simulation) 

O N=10,K=5(analytical approximation) 



Rate threshold -) (bits) 


Fig. 5. The comparison of the rate coverage probability obtained by Monte- 
Carlo simulations and the analytical approximation Gi. It is assumed that 
/3 = 4, G {3, 4,10}, and a: G {1, 2, 5}. 


in Fig. H] (b) in Section VI. As illustrated in Fig. |5] the 
analytical approximation tightly matches the simulated rate 
coverage probability over entire rate threshold region and for 
different numbers of BS antennas and users settings. 

One noticeable point in the rate coverage analysis is that 
since we do not assume any particular random network, 
the rate coverage is a function of distances between BSs 
and the tagged user. Prior work also characterized 

performance in a fixed network model, where performance 
was derived as a function of the distances between BSs and a 
user. 


IV. Ergodic Spectral Efficiency Analysis 

In this section, we derive an exact expression of the ergodic 
spectral efficiency in an integral form and provide a lower 
bound on the ergodic spectral efficiency in a closed form. As 
in the rate coverage analysis, the ergodic spectral efficiency is 
characterized as a function of distances between BSs and the 
tagged user. 


A. Exact Characterization 

Our exact expression is derived by using Lemma |2] which 
yields a general expression of the ergodic spectral efficiency 
in terms of moment generating functions. 


Lemma 2. Let xi,..., xat, i/i,..., um be arbitrary non-negative 
random variables. Then 


E 


AT 
2 -^n—l 


In 1 + 


E m , -I I 

m=l Vm + 1 / 

(^) (^) 


exp (—z) dz, 


(18) 


and M.x,y (z) = 


where Aiy (z) = E e 


Proof: See Lemma 1 in reference 
Leveraging Lemma |2l an exact expression of the ergodic 
spectral efficiency is obtained as in the following Theorem. 


Theorem 3. When the BSs form BS clusters by using the 
proposed method, the ergodic spectral efficiency of the tagged 
user is 


( lldoll'’ \ 
log2 (e) f°° V ^SNR/A ) 

L Jo ^ 

n 


K 


1 + z 


/JIM' 

I lldoll 


-P 


1 - 


/ ^ \ N—2K-\-l 


dz. 


(19) 


Proof: We start by characterizing the moment generating 
function of the desired channel gain. Since the desired channel 
gain is a Chi-squared random variable with 2 (V — 2K + 1) 
degrees of freedom from Lemma [T] its moment generating 
function is given by 


Ms {z) =E 


1 + z 


N-2K+1 


( 20 ) 


Similarly, the moment generating function of the interference 
is 



1 

l + z{||d,||/||do||)-^ 


K 


( 21 ) 


By applying Lemma |2l the Theorem is completed. ■ 

Given distances from the tagged user to each BS, the 
ergodic spectral efficiency of the tagged user is obtained by 
numerically calculating ( fTOl l. 


B. A Lower Bound Characterization 

We also derive a lower bound on the ergodic spectral 
efficiency with a closed form which is useful in obtaining 
insight on how system parameters change the ergodic spectral 
efficiency. To derive a lower bound, the following Lemma is 
claimed. 


Lemma 3. Lor any non-negative independent random vari¬ 
ables S and I, 


E 




> logs 


gE[lnS] \ 

E[/+1]J ■ 


( 22 ) 


Proof: 


E 




>Es 


E 


l 0 g 2 



(a) 

> E 




glnS Y 

E[/+l] 



gE[lnS] \ 

E[/+l] j 



(23) 


where (a) follows Jensen’s inequality and log 2 (l + 2) is 
a convex function of non-negative variable x, and (b) also 
follows Jensen’s inequality and log 2 (1 + e^) is a convex 
function of non-negative variable y. ■ 
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- N=4,K=1 (simulation) 

□ N=4,K=1(analytical lower bound) 

- N=3.K=1 (simulation) 

^ N=3,K=1(analytical lower bound) 

-N=4.K=2(simulation) 

^ N=4.K=2(analytical lower bound) 

- N=10,K=5(simulation) 

0 N=10,K=5(analytical lower bound) 

□ □ i 




* * * * ; 


// y* 

/a 0 0 0 ' 


/ y* 


50 60 70 80 90 100 110 120 

SNR(dB) 


Fig. 6. The comparison of the ergodic spectral efficiency obtained by Monte- 
Carlo simulations and the analytical lower bound m. It is assumed that 
/3 = 4, G {3, 4,10}, and a: G {1, 2, 5). 


Applying Lemma [3 a lower bound of the ergodic spectral 
efficiency is obtained in the following Theorem. 

Theorem 4. When the BSs form a BS coordination set 
according to the proposed cluster pattern, the ergodic spectral 
efficiency of the tagged user is lower bounded as 

i?^(SNR,||do||,P£,iV,iT,L,/3) 

exp {ijj {N — 2K + 1)) 


> f- log 2 


1 H- _p 

\ +l|dof/SNR/ 7 T^ 

(24) 


where V' {z) is the digamma function defined as T' (z) /T (z), 
where T (z) is the gamma function defined as 


/*00 

r(z)= / x^-^e-^dz. 

Jo 


Proof: First, we characterize E 


{2{N-2K + l)), 


In 


. Since 


hi 


E 


In 


= 'ijj{N -2K+ 1), 


(25) 
2 

(26) 


V. Ergodic Spectral Efficiency in Random 
Networks 

So far, we characterize the performance when the BSs 
are coordinated with the proposed cluster patterns in a fixed 
network model. In this section, we analyze the ergodic spectral 
efficiency when applying the proposed method in a random 
network. We skip characterizing the rate coverage probability 
of the proposed method in a random network since it is 
equivalent to a special case of (26) and (27) in ED. We 
consider a downlink cellular network where BSs equipped 
with N antennas are distributed according to a homogeneous 
PPP, 4> = {dijisN} with density A on the plane Eor 
easy of notation, we start the BS location index with 1 so 
that the fc-th closest BS from the origin is denoted as d^. 
Single antenna users are also distributed as a homogeneous 
PPP, <l)u = { Ui I z S N} with density Au, which is independent 
with <!>. We assume that Au ^ A so that there exist at 
least K users in each region with high enough probability. 
Under this assumption, K users are selected in each region. 
In this network model, the proposed method is applied in 
the same way with the case of a fixed network model, i.e., 
make BS clusters by using 2nd-order Voronoi region, draw a 
graph G (4)) by Delaunay triangulation, and design BS cluster 
patterns by using the corresponding edge-coloring. Under the 
assumption that the proposed method is used, we derive a 
lower bound on the ergodic spectral efficiency in a closed 
form. One should note that the randomness of the network 
model in this section is not a core part of the proposed method. 
Rather, it is a tool for analyzing the ergodic spectral efficiency 
performance. 

Without loss of generality, we focus on the typical user 
located at the origin, and assume that it is associated with the 
BS located at di. Under the premise that the typical user is 
in V 2 (di, d 2 ) G Vi, the typical user is served by a BS cluster 
{di, d 2 }. Similarly to the fixed network case, the typical user 
communicates with a BS pair {di,d 2 } by applying CBE, so 
that the intra-cluster interference from the BS at d 2 to the 
typical user is removed. Then, equivalently with (|9]), the SINR 
of the typical user is given by 


where (z) is the digamma function defined as E' (z) /E (z). 
Next, we compute E[J^]. Since the interference is a sum 
of independent Chi-squared random variables with different 
weights corresponding to the ratio of the distances, we have 

By applying Lemma Q, we complete the proof. ■ 

Now we verify the derived lower bound by comparing with 
the ergodic spectral efficiency obtained by using Monte-Carlo 
simulations. We use a network model described in Fig. |8] (b) 
in Section VI for the verification. It is assumed that L = 
Aec(G) = 4. The comparison result is in Fig. |6] In Fig. |6] 
we observe two separate regime depending on SNR, which 
are referred as the DoF regime and the saturation regime in 
0 - In both regimes, the lower bound is reasonably tight for 
different numbers of BS antennas and users settings. 


SINR|^ 


lldiir'' 

h{ 

2 

SdjG<I>\Z3(0.||d2||) dj 1 ^ 

(h^)"v^ 

^ -f 1 /snr/a: 


(28) 


where h{ is the modified channel coefficient from the BS at 
di to the typical user after applying CBF, is the channel 
coefficient vector from the BS at dj to the typical user, and 
is a beamforming matrix of the BS at dj, respectively. 
The entries of the channel coefficient vector follows Cfif (0,1). 
From (ED, the ergodic spectral efficiency is defined as 

(SNR, N, K, P) = ^E [log2 (l + SINR,,)] , (29) 

where the pre-log term 1/L is used because 1/L time- 
frequency resources are exploited to serve the typical user. 
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A. Lower Bound on the Ergodic Spectral Efficiency 


of the denominator of the SIR as 


Now we derive a lower bound on the ergodic spectral 
efficiency. To do this, we need the following Lemma. Lemma 
m provides probability density function (PDF) of distances of 
the closest BS and fc-th closest BS to the origin. 

Leituna 4. The joint probability density function (PDF) of 
||di|| and ||dfe|| is 


/||di||,||d*,|| {ri,rk) 

= [ H 

1 0 


k-2 


— Xirri 


if ri < Tk 
otherwise. 


(30) 


Proof: See Appendix C in ET\ . ■ 

Leveraging Lemma IH Theorem |5] is claimed. For analytical 
tractability, we assume the interference limited regime. 


E 


|dif E l|d* 

diG<l>\t5(0.||d2||) 




(a) 

= KE 


=KEr 


(t>) 

= KEr„r 
(c) 2KttX 


|dif 11 ^* 

diG<l>\t5(0.||d2||) 


E, 


■I'\S(0.r2) 


27rA 


E 


diG$\B( 0 ,r 2 ) 


l|di|| = ri, 

I|d2|| = r2 


^1-/3 dr 


nOO nr 2 

/ / 4(A7r)^e-^^’'=rf+V^"'^dridr2 

J r -}—0 J r ^ =0 


P 2 J 7ri —0 

SK 

' 52-4 


(35) 


Theorem 5. When BSs’ locations are distributed by a homo¬ 
geneous PPP and the proposed method is applied, the ergodic 
spectral efficiency of the typical user is lower bounded as 

> ;^log2 [l + -^^^^exp(Vi (iV-2iT+ !))], (31) 


where tp [z) is the digamma function, defined as F' {z) /F {z), 
where F {z) is the gamma function defined as 

pOO 

r{z)= / x^-^e-^dz. (32) 

JQ 

Proof: The ergodic spectral efficiency (l29l l is written by 


E 


i log2 (1 + SINR|,) 
^E [log 2 (1 + SINR|,)] 


(o) 1 

> -p l 0 g 2 



eE[in(|i.^ir)] 

A 



(33) 


where A = E ||dif Ed,e4.\t5(o.||d2||) l|d*ll 

(a) follows Lemma [3 and L = Aec(G)- Note that Aec(G) 

is determined as an algorithm parameter of the edge-cutting 

algorithm. Now, we obtain the expectation in the nominator 

and denominator of (l33l l. Since the interference is mitigated 

by CBF, as in the case of the deterministic square grid 
2 


network model, 
leads to 


hi 


follows (2 (N — 2K + 1)), which 



f / 


2 \ j 

E 

In f 

h[ 

)1 


f{N-2K + l), 


(34) 


where ip {■) is the digamma function. Next, we calculate the 
expectation of the interference part. We write the expectation 


where (a) follows ^ (2Ff) and it is indepen¬ 

dent to <i>, (b) follows the Campbell’s theorem, and (c) follows 
that the joint PDF in Lemma H] Plugging ( l34b and dTSl l into 
the proof is completed. ■ 

One should note that Theorem |5] considers the worst case 
for the out-of-cluster interference. This is because, in an 
irregular and asymmetric network model, e.g., a homogeneous 
PPP model, some BSs do not use the specific time-frequency 
resource exploited to serve the typical user. Then, the inter¬ 
ference from those BSs should not be counted in the out-of- 
cluster interference term. Analyzing this, however, is highly 
difficult because it depends on the correlation of the distributed 
points. It is not tractable in the known PPP framework. For 
this reason, we rather consider the worst case scenario, which 
assumes that the out of cluster interference comes from all the 
BSs in a network. 


VI. Performance Comparisons 
A. Simulation Setup 

We provide performance comparison results to demonstrate 
the superiority of the proposed method via system level simu¬ 
lations. As a simulation setup, we consider a 1400m x 1400m 
square network plane, that consists of 49 (= 7 x 7) square 
grids whose a side length is 200m. Within each square grid, 
one BS is randomly located. A range of a BS location is 
restricted by the perturbation square whose a side length is 
p, 0 < p < 200, and a center is same with the square 
grid. Given p, a BS can randomly located inside pm x pm 
perturbation square. For example, if the perturbation square 
has p = 0, each BS is regularly located on a center of each 
square grid. This is equivalent to a regular network. In other 
extreme case, if the perturbation parameter has p = 200, each 
BS can be located in any place in each square grid, which 
causes maximal randomness in a network. Fig. |7] describes 
the perturbation square and the corresponding BS location. 
This setting mimics an actual BS deployment scenario, where 
a BS is built to satisfy a certain guard region to other BSs, 
preventing that any two BSs are closely located each other. 
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200 m 



Fig. 7. A description of the perturbation square whose a side length is p. In 
the described case, a BS is randomly located inisde a region marked by grey 
shade. 



Fig. 8. Network snapshots for p = 100 in (a) and p = 200 in (b). It is 
assumed that A'perBS = 10. Each network plane is tessellated by 1st order 
Voronoi region of each BS. 


For the simulation, we consider two cases, p = 100 and 
p = 200. A case of p = 200 represents more irregularity of 
BSs locations than a case of p = 100. In these networks, 
ATperBS X 49 users are uniformly dropped. Network snap¬ 
shots for these two cases are illustrated in Fig. 0 For these 
snapshots, FfperBS = 10. Clearly, we observe that BSs are 
more irregularly located when p = 200. In the performance 
comparison, we only consider users inside 1000m x 1000m for 
wrap-around. The number of BS antennas and selected users 
are assumed to be = 3 and K = 1. The path-loss exponent 
is 4. 

Applying the proposed method in a given network, it is 
set that L = AecCG) = 4 in both cases of p = 100 
and p = 200 for the edge cutting algorithm. Recalling the 
edge cutting algorithm, we cut edges inversely proportional 
to the area of the corresponding 2nd-order Voronoi regions. 
However, measuring the area of 2nd-order Voronoi region 
is a challenging task. For this reason, instead of directly 
measuring each area, we uniformly drop dummy users and 
count the number of users associated with the certain two 
closest BSs pair. Denoting that the number of users whose the 
two closest BSs are {di,dj} as A(di,dj), we approximate 
A(V2(di,dj)) w A(di,dj). In the simulation, we drop 5000 
dummy users for the approximation. 

For the performance comparison, we consider two conven¬ 
tional coordination methods: dynamic clustering and static 
clustering. We explain as follows: 

1) Dynamic Clustering: In this method, selected users 
communicate with their two closest BSs. For comparison, 
we propose a simple user scheduling method by modifying 
FlashLinQ ITO) for dynamic clustering. Assuming a global 
scheduler, it assigns random orders to dropped users and the 


orders can be changed over time for user fairness. Then, 
according to the assigned order, a user finds the two closest 
BSs and create a BS cluster. If two users select the same BS, 
i.e., the BS selection conflict problem occurs, a user whose an 
order is late gives up communication in this particular time 
slot and wait for the next time slot. Since two users are able 
to communicate with one BS cluster (K = 1 for one BS and 
a BS cluster has two BSs), once a BS cluster is formed by 
a user, a global scheduler immediately finds the other user to 
communicate with the corresponding BS cluster. This process 
continues until the number of made BS clusters is maximal. 
Once all the BS clusters are made, selected users communicate 
with their BS clusters by using CBF, same with the proposed 
method. If there are remaining BSs that could not make a BS 
cluster, a global scheduler finds a user to communicate with 
them. No BS coordination is used for these users. 

2) Static Clustering: In this method, BS clusters are pre¬ 
determined irrespective of users’ conditions. A BS randomly 
selects one BS among its neighbor BSs, and makes a BS 
cluster with the selected one. Then, a user is independently 
selected by each BS and communicates with the predetermined 
BS cluster. Similarly to the proposed method and dynamic 
clustering, the intra-cluster interference is mitigated by using 
CBF. In the static clustering, there is non-zero probability that 
a user cannot communicate with its two closest BSs since BS 
clusters are made independently to users’ locations. 


B. Simulation Results 

First, we compare the edge users sum throughput depending 
on iTperBS- An edge user is defined as a user whose a distance 
ratio between the closest BS and the second closest BS is 
larger than 2/3, i.e., ||do|| / ||di || > 2/3. For obtaining the 
edge users sum throughput, we only pick the active edge users 
and sum their throughput. In the simulation, SNR is set to 
be lOOdB. Fig. |9] illustrates the edge users sum throughput. 
As observed in the figure, when the dropped users are dense 
enough, the performance of the proposed method achieves 
the performance of the dynamic clustering. Specifically, the 
proposed method have the same performance as the dynamic 
clustering when ATperBS > 20 in the p = 100 case, and 
when ATperBS > 40 in the p = 200 case. Comparing to 
the performance of the static clustering in the dense user 
environment, the edge users throughput is improved by 28% 
in the p = 100 case and 24% in the p = 200 case. The 
rationale behind this improvement obviously comes from that 
users communicate their two closest BSs in the proposed 
method, where the strongest interference is mitigated. In the 
static clustering, the strongest interference still can degrade the 
performance. In contrary to that, when the dropped users are 
not dense enough, the performance of the proposed method is 
even lower than that of the static clustering. This is because 
the proposed method uses predefined cooperative regions 
(2nd-order Voronoi regions) and time-frequency allocation. 
For intuition, let’s consider an example cooperative region 
1^2 (di, dj). For serving users in this region, a BS cluster 
{d„d j} is made and a specific time-frequency resource is 
allocated to this cluster. If there is no user in this region. 
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(a) (b) 

Fig. 9. The edge users sum throughput depending on iCperBS the p = 100 
case (a) and in the p = 200 case (b). It is assumed that N = 3, K = 1, and 
/3 = 4. The edge cutting is performed with Aec(G) = 4. 



Fig. 10. The rate coverage probability in the p = 100 case (a) and the 
p = 200 case (b). It is assumed that N = 3, K = 1, and /J = 4. The edge 
cutting is performed with Aec(G) = 4. 


however, the allocated time-frequency resource is wasted, and 
this degrades the performance. In the dynamic clustering, this 
resource waste does not occur since a BS cluster is flexibly 
made by a global scheduler, therefore there is always a user to 
communicate with the BS cluster. When the users are dense 
enough, all the predefined region has at least K users therefore 
no resource is wasted. Then, the proposed method provides the 
same performance as the dynamic clustering. This intuition 
also justifies the difference between the results of the p = 100 
and p = 200 cases. Since the BSs locations are more regular 
when p = 100, each 2nd-order Voronoi region has similar 
area. When p = 200, however, the BSs locations are highly 
irregular, therefore some 2nd-order Voronoi regions can be 
very small. Then, the probability that each region has at least 
K users is reduced. This is the reason why in the p = 200 
case, more users are required to achieve the same performance 
as the dynamic clustering compared to the p = 100 case. 

From this observation, we make some intuitive explanations 
about the effect of Aec(G) on the network performance. 
Considering Aec(G') = 1 as an extreme case, no time- 
frequency resource are wasted. A large user density, however, 
would be required to achieve the same level of the coordination 
benefit with that of dynamic clustering. This is mainly because 
many 2nd-order Voronoi regions are cut, therefore users exist 
in the remaining 2nd-order Voronoi region with only small 
probability. In another extreme case, if the edge-cutting algo¬ 
rithm is not applied, i.e., Aec(G) = A(G), the performance 
of the proposed method is not affected by the user density 
since all regions are covered by the proposed method, while 
the network performance degradation is unavoidable due to 
the resource waste. In summary, less user density is required 
as Aec (G) increases, while more resource waste is expected. 

We can incorporate the channel estimation overhead into the 
sum throughput performance. From im, the pilot overhead is 

Lp(2,A,SINR) 

Of , (3u) 

where Lp{2, N, SINK) is the number of symbols for pilots and 
Lb is a given fading coherence time in symbols. The number 
of symbols for pilots is Lp(2, N, SINK) = 2r]N, and rj is 


rj = max 


1 , 


1 

SINK 


(-L_ 

Vmmse 



(37) 


where MMSE is the minimum mean square error of the pilot 
signals. Incorporating the pilot overhead into the throughput. 


the performance is degraded by (1 — a). Since all the con¬ 
sidered coordination methods are pair-wise BS coordination 
where a BS cluster consists of two BSs, the amount of channel 
estimation overhead is same. 

Now we compare the rate coverage probability. For obtain¬ 
ing the rate coverage probability, we calculate the probability 
that the selected user’s throughput is larger than rate threshold 
7. It is assumed that SNR = lOOdB and iTperBS = 40. 
Fig. Ho] illustrates the rate coverage probability. As observed in 
this figure, the proposed method provides the same or similar 
coverage probability with that of the dynamic clustering at 
low rate threshold, while the coverage probability degrades to 
the similar level of the static clustering at high rate threshold. 
Focusing on 7 = Ibit, the rate outage probability is reduced 
by 50% comparing to the static clustering in the p = 100 
case, and it is reduced by 15% in the p = 200 case. 
Since the rate coverage performance at low rate threshold is 
determined by edge users, the simulation results show that 
the proposed method is able to provide the same or similar 
level of coordination benefit with the dynamic clustering to 
edge users. At high rate threshold, however, the rate coverage 
performance is degraded. The reason of this observation is on 
the edge cutting algorithm. As a result of the edge cutting 
algorithm, we cut some regions and do not serve users in the 
cut regions even if users’ locations in the cut region are very 
close to their associated BS. Then, we miss a chance to have 
a high rate, which degrades the rate coverage at high rate 
threshold. 


VII. Conclusion 

In this paper, we proposed a method for designing BS 
clusters and BS cluster patterns in irregular BS topologies. 
The core idea of the proposed method is to form BS clusters 
according to the 2nd-order Voronoi regions, and use edge¬ 
coloring for the graph induced by Delaunay triangulation to 
assign time-frequency resources (colors) to each BS cluster. 
By using the proposed method, users in the uncut regions 
communicate with their two closest BS while avoiding the BS 
selection conflict problem. In a fixed network with irregular 
BSs’ locations, we provided analytical expressions of the rate 
coverage probability and the ergodic spectral efficiency in 
terms of the number of antennas per BS, the number of users, 
distance, SNR, and the pathloss exponent. In a random net¬ 
work modeled by a homogeneous PPP, we also derived a lower 
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bound of the ergodic spectral efficiency as a function of the 
pathloss exponent and the number of associated users. Through 
system level simulations, we showed that the proposed method 
achieves the same edge users sum throughput performance 
with dynamic clustering when the users are dense enough. 
In the rate coverage probability performance, the proposed 
method provides the same or similar performance with that 
of dynamic clustering at low rate threshold. Future work 
should consider the application of more advanced coordination 
strategies going beyond coordinated beamforming. 
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